DNA-based geographic typing of the gourmet mushroom Tricholoma matsutake traded in China

Metals Transfer in Mushroom Tricholoma matsutake from Regional High Geochemical Background Areas: Environmental Influences and Human Health Risk

Wild-grown edible mushrooms are important in world diets and are also efficient metal accumulators. Yunnan, Southwest China, is the main producing region, with typically high levels of geochemical metals. The environmental factors, bioaccumulation, distribution and human health risks of metals were examined in paired soil and Tricholoma matsutake (n = 54). T. matsutake grows on acidified soils (pH = 3.95-6.56), and metals show a strong heterogeneity, with Fe, Mn, Zn and Cu in the ranges of 16-201, 0.046-8.58 g kg-1, and 22.6-215, 3.7-155 mg kg-1. High soil Fe content led to great accumulation in T. matsutake (0.24-18.8 g kg-1). However, though the soil Mn content was higher than that of Zn and Cu, their concentrations in T. matsutake were comparable (21.1-487 vs. 38.7-329 and 24.9-217 mg kg-1). This suggested that T. matsutake prefers to accumulate Zn and Cu compared to Mn, and this is supported by the bioaccumulation factors (BAFs = 0.32-17.1 vs. 0.006-1.69). Fe was mainly stored in stipes, while Mn, Zn and Cu were stored in caps, and the translocation factors (TFs) were 0.58 vs. 1.28-1.94. Therefore, stipe Fe showed the highest health risk index (HRI) at 1.28-26.9, followed by cap Cu (1.01-2.33), while 98-100% of the Mn and Zn were risk-free. The higher concentration and greater risk of Fe was attributed to the significant effect of soil Fe content (R = 0.34) and soil pH (R = -0.57). This study suggested that Fe, as an essential mineral, may exert toxic effects via the consumption of T. matsutake from high geochemical background areas.

Dynamic evolution of eukaryotic mitochondrial and nuclear genomes: a case study in the gourmet pine mushroom Tricholoma matsutake

Fungi, as eukaryotic organisms, contain two genomes, the mitochondrial genome and the nuclear genome, in their cells. How the two genomes evolve and correlate to each other is debated. Herein, taking the gourmet pine mushroom Tricholoma matsutake as an example, we performed comparative mitogenomic analysis using samples collected from diverse locations and compared the evolution of the two genomes. The T. matsutake mitogenome encodes 49 genes and is rich of repetitive and non-coding DNAs. Six genes were invaded by up to 11 group I introns, with one cox1 intron cox1P372 showing presence/absence dynamics among different samples. Bioinformatic analyses suggested limited or no evidence of mitochondrial heteroplasmy. Interestingly, hundreds of mitochondrial DNA fragments were found in the nuclear genome, with several larger than 500 nt confirmed by PCR assays and read count comparisons, indicating clear evidence of transfer of mitochondrial DNA into the nuclear genome. Nuclear DNA of T. matsutake showed a higher mutation rate than mitochondrial DNA. Furthermore, we found evidence of incongruence between phylogenetic trees derived from mitogenome and nuclear DNA sequences. Together, our results reveal the dynamic genome evolution of the gourmet pine mushroom.

Exploring the Species Diversity of Edible Mushrooms in Yunnan, Southwestern China, by DNA Barcoding

Yunnan Province, China, is famous for its abundant wild edible mushroom diversity and a rich source of the world's wild mushroom trade markets. However, much remains unknown about the diversity of edible mushrooms, including the number of wild edible mushroom species and their distributions. In this study, we collected and analyzed 3585 mushroom samples from wild mushroom markets in 35 counties across Yunnan Province from 2010 to 2019. Among these samples, we successfully obtained the DNA barcode sequences from 2198 samples. Sequence comparisons revealed that these 2198 samples likely belonged to 159 known species in 56 different genera, 31 families, 11 orders, 2 classes, and 2 phyla. Significantly, 51.13% of these samples had sequence similarities to known species at lower than 97%, likely representing new taxa. Further phylogenetic analyses on several common mushroom groups including 1536 internal transcribed spacer (ITS) sequences suggested the existence of 20 new (cryptic) species in these groups. The extensive new and cryptic species diversity in wild mushroom markets in Yunnan calls for greater attention for the conservation and utilization of these resources. Our results on both the distinct barcode sequences and the distributions of these sequences should facilitate new mushroom species discovery and forensic authentication of high-valued mushrooms and contribute to the scientific inventory for the management of wild mushroom markets.

Comparative analysis of transcriptomes revealed the molecular mechanism of development of Tricholoma matsutake at different stages of fruiting bodies

The purpose of the study is to investigate the molecular mechanisms of development of Tricholoma matsutake fruiting body at the primordial stage (TM-1), the intermediate stage (TM-2) and the mature stage (TM-3) using RNA-Seq sequencing technology. The analysis of gene expression level revealed that the Spn2 and Eef1a1 gene were the key genes in the primordial stage of T. matsutake by regulating cytokinesis, protein synthesis, and cell growth. And the Ubc, Atp6, Cytb, and Pth2 gene were the key genes in the mature stage of T. matsutake by regulating energy metabolism and protein synthesis. Differential expression genes (DEGs) analysis results showed that Cdc28, Rad53, Dun1, Pho85 and Pho81 were the key DEGs regulating cell cycle genes of T. matsutake from primordial stage to intermediate stage. And APC, Cyr1, Cdc45, Spo11 and Rec8 genes were the key DEGs for the meiosis and sporogenesis of T. matsutake from the intermediate stage to the mature stage.

Prediction of the potential geographic distribution of the ectomycorrhizal mushroom Tricholoma matsutake under multiple climate change scenarios

Effective conservation and utilization strategies for natural biological resources require a clear understanding of the geographic distribution of the target species. Tricholoma matsutake is an ectomycorrhizal (ECM) mushroom with high ecological and economic value. In this study, the potential geographic distribution of T. matsutake under current conditions in China was simulated using MaxEnt software based on species presence data and 24 environmental variables. The future distributions of T. matsutake in the 2050s and 2070s were also projected under the RCP 8.5, RCP 6, RCP 4.5 and RCP 2.6 climate change emission scenarios described in the Special Report on Emissions Scenarios (SRES) by the Intergovernmental Panel on Climate Change (IPCC). The areas of marginally suitable, suitable and highly suitable habitats for T. matsutake in China were approximately 0.22 × 10⁶ km², 0.14 × 10⁶ km², and 0.11 × 10⁶ km², respectively. The model simulations indicated that the area of marginally suitable habitats would undergo a relatively small change under all four climate change scenarios; however, suitable habitats would significantly decrease, and highly suitable habitat would nearly disappear. Our results will be influential in the future ecological conservation and management of T. matsutake and can be used as a reference for studies on other ectomycorrhizal mushroom species.

A qPCR assay that specifically quantifies Tricholoma matsutake biomass in natural soil

Tricholoma matsutake is an ectomycorrhizal basidiomycete that produces prized, yet uncultivable, "matsutake" mushrooms along densely developed mycelia, called "shiro," in the rhizosphere of coniferous forests. Pinus densiflora is a major host of this fungus in Japan. Measuring T. matsutake biomass in soil allows us to determine the kinetics of fungal growth before and after fruiting, which is useful for analyzing the conditions of the shiro and its surrounding mycorrhizosphere, predicting fruiting timing, and managing forests to obtain better crop yields. Here, we document a novel method to quantify T. matsutake mycelia in soil by quantifying a single-copy DNA element that is uniquely conserved within T. matsutake but is absent from other fungal species, including close relatives and a wide range of ectomycorrhizal associates of P. densiflora. The targeted DNA region was amplified quantitatively in cultured mycelia that were mixed with other fungal species and soil, as well as in an in vitro co-culture system with P. densiflora seedlings. Using this method, we quantified T. matsutake mycelia not only from shiro in natural environments but also from the surrounding soil in which T. matsutake mycelia could not be observed by visual examination or distinguished by other means. It was demonstrated that the core of the shiro and its underlying area in the B horizon are predominantly composed of fungal mycelia. The fungal mass in the A or A₀ horizon was much lower, although many white mycelia were observed at the A horizon. Additionally, the rhizospheric fungal biomass peaked during the fruiting season.