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.

Identification of the Pol Gene as a Species-Specific Diagnostic Marker for Qualitative and Quantitative PCR Detection of Tricholoma matsutake.

Tricholoma matsutake is a rare, precious, and wild edible fungus that could not be cultivated artificially until now. This situation has given way to the introduction of fake T. matsutake commodities to the mushroom market. Among the methods used to detect food adulteration, amplification of species-specific diagnostic marker is particularly important and accurate. In this study, the Pol gene is reported as a species-specific diagnostic marker to identify three T. matsutake varieties and 10 other types of edible mushrooms through qualitative and quantitative PCR. The PCR results did not reveal variations in the amplified region, and the detection limits of qualitative and quantitative PCR were found to be 8 ng and 32 pg, respectively. Southern blot showed that the Pol gene exists as a single copy in the T. matsutake genome. The method that produced the purest DNA of T. matsutake in this study was also determined, and the high-concentration salt precipitation method was confirmed to be the most suitable among the methods tested. The assay proposed in this work is applicable not only to the detection of raw materials but also to the examination of processed products containing T. matsutake.

Heavy-ion beam mutagenesis of the ectomycorrhizal agaricomycete Tricholoma matsutake that produces the prized mushroom "matsutake" in conifer forests.

Tricholoma matsutake is an ectomycorrhizal agaricomycete that produces the prized mushroom "matsutake" in Pinaceae forests. Currently, there are no available cultivars or cultivation methods that produce fruiting bodies. Heavy-ion beams, which induce mutations through double-stranded DNA breaks, have been used widely for plant breeding. In the present study, we examined whether heavy-ion beams could be useful in isolating T. matsutake mutants. An argon-ion beam gave a suitable lethality curve in relation to irradiation doses, accelerating killing at 100-150 Gy. Argon-ion beam irradiation of the agar plate cultures yielded several transient mutants whose colony morphologies differed from that of the wild-type strain at the first screening, but which did not persist following culture transfer. It also generated a mutant whose phenotype remained stable after repeated culture transfers. The stable pleiotropic mutant not only exhibited a different colony morphology to the wild type, but also showed increased degradation of dye-linked water-insoluble amylose and cellulose substrates. Thus, heavy-ion beams may be useful for isolating mutants of T. matsutake, although precautions may be required to maintain the mutants, without phenotypic reversion, during repetitive culture of their mycelia.

North American matsutake: names clarified and a new species described.

Tricholoma matsutake, known widely as "matsutake," has great commercial and cultural significance in Japan. Because Japanese production is insufficient to meet the high domestic demand, morphologically similar mushrooms, thought by many to belong to T. magnivelare, are imported from western North America. However, molecular data produced since the early 2000s have indicated that more than one species of matsutake occur in North America and this raises the question of correct naming for the different species. To address this question, we assessed the phylogenetic diversity within North American matsutake based on nuc rDNA ITS1-5.8S-ITS2 (internal transcribed spacer [ITS] barcode) sequences, including newly obtained sequences from the type collections for Agaricus ponderosus and Armillaria arenicola, and morphological characters. Our results agree with earlier indications that three matsutake species occur in North America and allow us to clarify the correct application of names-T. magnivelare from the eastern USA and Canada, T. murrillianum from the western USA and Canada, and T. mesoamericanum from Mexico, newly described here. The existence of the three North American species is further supported by the results of evolutionary divergence analysis, geographical distributions, and morphological characters.

Dehydrogenase genes in the ectomycorrhizal fungus Tricholoma vaccinum: A role for Ald1 in mycorrhizal symbiosis.

Ectomycorrhizal symbiosis is important for forest ecosystem functioning with tree-fungal cooperation increasing performance and countering stress conditions. Aldehyde dehydrogenases (ALDHs) are key enzymes for detoxification and thus may play a role in stress response of the symbiotic association. With this focus, eight dehydrogenases, Ald1 through Ald7 and TyrA, of the ectomycorrhizal basidiomycete Tricholoma vaccinum were characterized and phylogenetically investigated. Functional analysis was performed through differential expression analysis by feeding different, environmentally important substances. A strong effect of indole-3-acetic acid (IAA) was identified, linking mycorrhiza formation and auxin signaling between the symbiosis partners. We investigated ald1 overexpressing strains for performance in mycorrhiza with the host tree spruce (Picea abies) and observed an increased width of the apoplast, accommodating the Hartig' net hyphae of the T. vaccinum over-expressing transformants. The results support a role for Ald1 in ectomycorrhiza formation and underline functional differentiation within fungal aldehyde dehydrogenases in the family 1 of ALDHs.

Genetic variability and bottleneck detection of four Tricholoma matsutake populations from northeastern and southwestern China.

The excessive commercial collection of matsutake mushrooms can lead to extreme reduction of population size, which may cause genetic bottleneck and decrease genetic diversity of Tricholoma matsutake. Here, six polymorphic microsatellite loci markers were used to examine the genetic diversity of four natural T. matsutake populations from two main producing regions of China. The minimum combinations of four loci were able to discriminate total 86 sampled individuals with distinctive multilocus genotypes. Our analysis of molecular variance (AMOVA) revealed that about 80% and 20% of the overall genetic variation were respectively partitioned within and among populations. The principal-coordinate analyses (PCA) distinguished the four tested populations into three genetic clusters, each of which was correlated with respective endemic host plants on a geographical basis. The AMOVA, PCA and pairwise population FST estimates consistently displayed the same genetic divergence patterns and spatial structure of T. matsutake mediated by host plants in China. The significant heterozygosity excesses demonstrated that a recent genetic bottleneck occurred in each population tested. The complementary M-ratio test indicated past genetic bottleneck events over longer periods. Only four individuals were identified as putative first generation migrants within northeastern China, which implies restricted interpopulation gene flow in T. matsutake. We discuss that the significant genetic differentiation among populations of T. matsutake is most likely a function of host adaptation, host specificity, genetic bottleneck, limited dispersal and habitat fragmentation.

Biosynthesis and Secretion of Indole-3-Acetic Acid and Its Morphological Effects on Tricholoma vaccinum-Spruce Ectomycorrhiza.

Fungus-derived indole-3-acetic acid (IAA), which is involved in development of ectomycorrhiza, affects both partners, i.e., the tree and the fungus. The biosynthesis pathway, excretion from fungal hyphae, the induction of branching in fungal cultures, and enhanced Hartig net formation in mycorrhiza were shown. Gene expression studies, incorporation of labeled compounds into IAA, heterologous expression of a transporter, and bioinformatics were applied to study the effect of IAA on fungal morphogenesis and on ectomycorrhiza. Tricholoma vaccinum produces IAA from tryptophan via indole-3-pyruvate, with the last step of this biosynthetic pathway being catalyzed by an aldehyde dehydrogenase. The gene ald1 was found to be highly expressed in ectomycorrhiza and induced by indole-3-acetaldehyde. The export of IAA from fungal cells is supported by the multidrug and toxic extrusion (MATE) transporter Mte1 found in T. vaccinum. The addition of IAA and its precursors induced elongated cells and hyphal ramification of mycorrhizal fungi; in contrast, in saprobic fungi such as Schizophyllum commune, IAA did not induce morphogenetic changes. Mycorrhiza responded by increasing its Hartig net formation. The IAA of fungal origin acts as a diffusible signal, influencing root colonization and increasing Hartig net formation in ectomycorrhiza.

Mobile DNA distributions refine the phylogeny of "matsutake" mushrooms, Tricholoma sect. Caligata.

"Matsutake" mushrooms are formed by several species of Tricholoma sect. Caligata distributed across the northern hemisphere. A phylogenetic analysis of matsutake based on virtually neutral mutations in DNA sequences resolved robust relationships among Tricholoma anatolicum, Tricholoma bakamatsutake, Tricholoma magnivelare, Tricholoma matsutake, and Tricholoma sp. from Mexico (=Tricholoma sp. Mex). However, relationships among these matsutake and other species, such as Tricholoma caligatum and Tricholoma fulvocastaneum, were ambiguous. We, therefore, analyzed genomic copy numbers of σ marY1 , marY1, and marY2N retrotransposons by comparing them with the single-copy mobile DNA megB1 using real-time polymerase chain reaction (PCR) to clarify matsutake phylogeny. We also examined types of megB1-associated domains, composed of a number of poly (A) and poly (T) reminiscent of RNA-derived DNA elements among these species. Both datasets resolved two distinct groups, one composed of T. bakamatsutake, T. fulvocastaneum, and T. caligatum that could have diverged earlier and the other comprising T. magnivelare, Tricholoma sp. Mex, T. anatolicum, and T. matsutake that could have evolved later. In the first group, T. caligatum was the closest to the second group, followed by T. fulvocastaneum and T. bakamatsutake. Within the second group, T. magnivelare was clearly differentiated from the other species. The data suggest that matsutake underwent substantial evolution between the first group, mostly composed of Fagaceae symbionts, and the second group, comprised only of Pinaceae symbionts, but diverged little within each groups. Mobile DNA markers could be useful in resolving difficult phylogenies due to, for example, closely spaced speciation events.

Near-complete de novo assembly of Tricholoma bakamatsutake chromosomes revealed the structural divergence and differentiation of Tricholoma genomes.

Tricholoma bakamatsutake, which is an edible ectomycorrhizal fungus associated with Fagaceae trees, may have diverged before the other species in Tricholoma section Caligata. We generated a highly contiguous whole-genome sequence for T. bakamatsutake SF-Tf05 isolated in an Oak (Quercus salicina) forest in Japan. The assembly of high-fidelity long reads, with a median read length of 12.3 kb, resulted in 13 chromosome-sized contigs comprising 142,068,211 bases with an average guanine and cytosine (GC) content of 43.94%. The 13 chromosomes were predicted to encode 11,060 genes. A contig (122,566 bases) presumably containing the whole circular mitochondrial genome was also recovered. The chromosome-wide comparison of T. bakamatsutake and Tricholoma matsutake (TMA_r1.0) indicated that the basic number of chromosomes (13) was conserved, but the structures of the corresponding chromosomes diverged, with multiple inversions and translocations. Gene conservation and cluster analyses revealed at least 3 phylogenetic clades in Tricholoma section Caligata. Specifically, all T. bakamatsutake strains belonged to the "bakamatsutake" clade, which is most proximal to the "caligatum" clade consisting of Tricholoma caligatum and Tricholoma fulvocastaneum. The constructed highly contiguous nearly telomere-to-telomere genome sequence of a T. bakamatsutake isolate will serve as a fundamental resource for future research on the evolution and differentiation of Tricholoma species.

Intercontinental divergence in the Populus-associated ectomycorrhizal fungus, Tricholoma populinum.

The ectomycorrhizal fungus Tricholoma populinum is host-specific with Populus species. T. populinum has wind-dispersed progagules and may be capable of long-distance dispersal. In this study, we tested the hypothesis of a panmictic population between Scandinavia and North America. DNA sequences from five nuclear loci were used to assess phylogeographic structure and nucleotide divergence between continents. Tricholoma populinum was composed of Scandinavian and North American lineages with complete absence of shared haplotypes and only one shared nucleotide mutation. Divergence of these lineages was estimated at approx. 1.7-1.0 million yr ago (Ma), which occurred after the estimated divergence of host species Populus tremula and Populus balsamifera/Populus trichocarpa at 5 Ma. Phylogeographic structure was not observed within Scandinavian or North American lineages of T. populinum. Intercontinental divergence appears to have resulted from either allopatric isolation; a recent, rare long-distance dispersal founding event followed by genetic drift; or the response in an obligate mycorrhizal fungus with a narrow host range to contractions and expansion of host distribution during glacial and interglacial episodes within continents. Understanding present genetic variation in populations is important for predicting how obligate symbiotic fungi will adapt to present and future changing climatic conditions.

Phylogenetic relationship and species delimitation of matsutake and allied species based on multilocus phylogeny and haplotype analyses.

Tricholoma matsutake (S. Ito & S. Imai) Singer and its allied species are referred to as matsutake worldwide and are the most economically important edible mushrooms in Japan. They are widely distributed in the northern hemisphere and established an ectomycorrhizal relationship with conifer and broadleaf trees. To clarify relationships among T. matsutake and its allies, and to delimit phylogenetic species, we analyzed multilocus datasets (ITS, megB1, tef, gpd) with samples that were correctly identified based on morphological characteristics. Phylogenetic analyses clearly identified four major groups: matsutake, T. bakamatsutake, T. fulvocastaneum and T. caligatum; the latter three species were outside the matsutake group. The haplotype analyses and median-joining haplotype network analyses showed that the matsutake group included four closely related but clearly distinct taxa (T. matsutake, T. anatolicum, Tricholoma sp. from Mexico and T. magnivelare) from different geographical regions; these were considered to be distinct phylogenetic species.

Modulation of ethanol stress tolerance by aldehyde dehydrogenase in the mycorrhizal fungus Tricholoma vaccinum.

We report the first mycorrhizal fungal aldehyde dehydrogenase gene, ald1, which was isolated from the basidiomycete Tricholoma vaccinum. The gene, encoding a protein Ald1 of 502 amino acids, is up-regulated in ectomycorrhiza. Phylogenetic analyses using 53 specific fungal aldehyde dehydrogenases from all major phyla in the kingdom of fungi including Ald1 and two partial sequences of T. vaccinum were performed to get an insight in the evolution of the aldehyde dehydrogenase family. By using competitive and real-time RT-PCR, ald1 is up-regulated in response to alcohol and aldehyde-related stress. Furthermore, heterologous expression of ald1 in Escherichia coli and subsequent in vitro enzyme activity assay demonstrated the oxidation of propionaldehyde and butyraldehyde with different kinetics using either NAD(+) or NADP(+) as cofactors. In addition, overexpression of ald1 in T. vaccinum after Agrobacterium tumefaciens-mediated transformation increased ethanol stress tolerance. These results demonstrate the ability of Ald1 to circumvent ethanol stress, a critical function in mycorrhizal habitats.

Effect of heavy metal-solubilizing microorganisms on zinc and cadmium extractions from heavy metal contaminated soil with Tricholoma lobynsis.

The macrofungus, Tricholoma lobynsis, was chosen to remedy Zn-Cd-Pb contaminated soil. To enhance its metal-extracting efficiency, two heavy metal resistant microbes M6 and K1 were applied owing to their excellent abilities to solubilize heavy metal salts. The two isolated microbial strains could also produce indole acetic acid (IAA), siderophore and solubilize inorganic phosphate, but neither of them showed 1-aminocyclopropane-1-carboxylate deaminase activity. The strains M6 and K1 were identified as Serratia marcescens and Rhodotorula mucilaginosa based on 16S rDNA and ITS sequence analysis respectively. Pot experiment showed that spraying to T. lobynsis-inoculated soil with M6 and K1 respectively could increase total Cd accumulations of this mushroom by 216 and 61%, and Zn by 153 and 49% compared to the uninoculated control. Pb accumulation however, was too low (<1 mg kg(-1)) to be determined. The results illustrated that special microbes and macrofungi can work together to remedy polluted soil as plant and plant growth promoting microbes do, probably because of excellent metal-accumulating abilities of macrofungi and IAA-siderophore production, phosphate solubilization abilities of the assisted-microbes. This kind of macrofungi-microbe interaction can be developed into a novel bioremediation strategy.

Quantitative transcriptomic and metabolomic analyses reveal the changes in Tricholoma matsutake fruiting bodies during cold storage.

A combination of transcriptomic and metabolomic analyses was performed to systematically understand the metabolic changes in Tricholoma matsutake fruiting bodies during cold storage. In total, 800 metabolites were identified and 19,964 annotated unigenes were quantified. The unigenes related to the catabolism of proteins, carbohydrates, and lipids were mainly upregulated during cold storage, but the related primary metabolites were not accumulated, which indicated complete degradation and loss of nutrients. Concurrently, the synthesis and metabolism of the main components of the cell wall, chitin and ß-1,3-glucan, were regulated, indicating the dynamic remodeling of the T. matsutake cell wall structure. Additionally, indole-3-acetic acid and components of its synthesis pathway were found in T. matsutake, indicating their potential role as a communicator between T. matsutake and its symbiotic plants. The results provide new information to improve the understanding of the metabolic mechanism of T. matsutake fruit bodies during postharvest cold storage.

The ectomycorrhizal association of Tricholoma matsutake and two allied species, T. bakamatsutake and T. fulvocastaneum, with native hosts in subtropical China.

Accurate identification of edible ectomycorrhizal (ECM) mushrooms and their host trees in nature is key to commercial production for consumption. For the first time we describe the ectomycorrhizal association of the three most common species of edible matsutake mushrooms with their native host trees in Yunnan Province in China. We collected ECM samples from three different localities in subtropical forests known to be highly productive areas of Tricholoma. Additionally, we collected basidiomata of Tricholoma from the field and markets in Yunnan. ECM samples were analyzed using morphological and molecular methods. We conducted phylogenetic analyses of nuc rDNA internal transcribed spacer region ITS1-5.8S-ITS2 (ITS) and analyzed the intergenic spacer of cpDNA psbA-trnH to identify basidiomata and plant hosts, respectively. Three species of Tricholoma were identified: T. bakamatsutake, T. fulvocastaneum, and T. matsutake. Four ECM associations in the study area were detected: Tricholoma bakamatsutake + Castanopsis tibetana, T. fulvocastaneum + C. tibetana, T. fulvocastaneum + Pinus yunnanensis, and T. matsutake + P. yunnanensis. Detailed descriptions and illustrations of the ECM associations are presented.

Tricholoma alpinum, sp. nov., under five-needle pines in alpine and subalpine zones in Japan.

Alpine regions in Japan are characterized by the Siberian dwarf pine, Pinus pumila. Although these regions are conserved as national parks due to their unique biome, few reports of the fungal biota are available. We examined mushroom fungi under such vegetation and found a unique yellowish-capped Tricholoma species. Multilocus molecular phylogenetic analyses and morphological observation of specimens showed that the Tricholoma species is very similar to T. fumosoluteum known in North America. The Japanese yellowish-capped T. aff. fumosoluteum had larger basidiospores and basidia and shorter pileipellis hyphae compared with T. fumosoluteum. Therefore, we describe the Japanese entity as a new species, T. alpinum.

Identification of upregulated genes in Tricholoma matsutake mycorrhiza.

Many plant roots associate with fungi to form mycorrhizae; tree roots especially associate with ectomycorrhizal fungi, such as Tricholoma species. Tricholoma matsutake is an economically important fungus in Asian countries and usually inhabits forests primarily composed of Pinus densiflora (Japanese red pine). In this study, to understand the mycorrhizal association between T. matsutake and P. densiflora, genes specifically expressed in mycorrhiza compared with those expressed in mycelia and fruiting bodies were identified by RNA-seq. This revealed that genes for chromatin, proteasomes, signal transduction, pheromones, cell surface receptors, cytoskeleton, RNA processing and transporters from T. matsutake were highly expressed in mycorrhiza. It also identified 35 mycorrhiza-induced small secreted proteins (MiSSPs) that were highly expressed in mycorrhiza. Meanwhile, genes for proteases, defence-related proteins, cell-wall degradation, signal transduction, pinene synthesis, plant hormones and transporters from P. densiflora were highly expressed in mycorrhiza. These genes may be involved in mycorrhizal formation and maintenance. A MiSSP, 1460819, was highly expressed in mycorrhiza, and this expression was maintained for 24 months. These results provide insight into the mycorrhizal association between T. matsutake and P. densiflora.

Phylogenetic species delimitation in ectomycorrhizal fungi and implications for barcoding: the case of the Tricholoma scalpturatum complex (Basidiomycota).

Population studies have revealed that the fungal ectomycorrhizal morphospecies Tricholoma scalpturatum consists of at least two genetically distinct groups that occur sympatrically in several geographical areas. This discovery prompted us to examine species boundaries and relationships between members formerly assigned to T. scalpturatum and allied taxa using phylogenetic analyses. Sequence data were obtained from three nuclear DNA regions [internal transcribed spacer (ITS), gpd and tef], from 101 carpophores collected over a large geographical range in Western Europe, and some reference sequences from public databases. The ITS was also tested for its applicability as DNA barcode for species delimitation. Four highly supported phylogenetic clades were detected. The two previously detected genetic groups of T. scalpturatum were assigned to the phylospecies Tricholoma argyraceum and T. scalpturatum. The two remaining clades were referred to as Tricholoma cingulatum and Tricholoma inocybeoides. Unexpectedly, T. cingulatum showed an accelerated rate of evolution that we attributed to narrow host specialization. This study also reveals recombinant ITS sequences in T. inocybeoides, suggesting a hybrid origin. The ITS was a useful tool for the determination of species boundaries: the mean value of intraspecific genetic distances in the entire ITS region (including 5.8S rDNA) was <0.2%, whereas interspecific divergence estimates ranged from 1.78% to 4.22%. Apart from giving insights into the evolution of the T. scalpturatum complex, this study contributes to the establishment of a library of taxonomically verified voucher specimens, an a posteriori correlation between phenotype and genotype, and DNA barcoding of ectomycorrhizal fungi.

Biosynthesis of (R)-(-)-1-Octen-3-ol in Recombinant Saccharomyces cerevisiae with Lipoxygenase-1 and Hydroperoxide Lyase Genes from Tricholoma matsutake.

Tricholoma matsutake is an ectomycorrhizal fungus, related with the host of Pinus densiflora. Most of studies on T. matsutake have focused on mycelial growth, genes and genomics, phylogenetics, symbiosis, and immune activity of this strain. T. matsutake is known for its unique fragrance in Eastern Asia. The most major component of its scent is (R)-(-)-1-octen-3-ol and is biosynthesized from the substrate linoleic acid by the sequential reaction of lipoxygenase and peroxide lyase. Here, we report for the first time the biosynthesis of (R)-(-)- 1-octen-3-ol of T. matsutake using the yeast Saccharomyces cerevisiae as a host. In this study, cDNA genes correlated with these reactions were cloned from T. matsutake, and expression studies of theses genes were carried out in the yeast Saccharomyces cerevisiae. The product of these genes expression study was carried out with Western blotting. The biosynthesis of (R)-(-)- 1-octen-3-ol of T. matsutake in recombinant Saccharomyces cerevisiae was subsequently identified with GC-MS chromatography analysis. The biosynthesis of (R)-(-)-1-octen-3-ol with S. cerevisiae represents a significant step forward.

Unusual genome expansion and transcription suppression in ectomycorrhizal Tricholoma matsutake by insertions of transposable elements.

Genome sequencing of Tricholoma matsutake revealed its unusually large size as 189.0 Mbp, which is a consequence of extraordinarily high transposable element (TE) content. We identified that 702 genes were surrounded by TEs, and 83.2% of these genes were not transcribed at any developmental stage. This observation indicated that the insertion of TEs alters the transcription of the genes neighboring these TEs. Repeat-induced point mutation, such as C to T hypermutation with a bias over "CpG" dinucleotides, was also recognized in this genome, representing a typical defense mechanism against TEs during evolution. Many transcription factor genes were activated in both the primordia and fruiting body stages, which indicates that many regulatory processes are shared during the developmental stages. Small secreted protein genes (<300 aa) were dominantly transcribed in the hyphae, where symbiotic interactions occur with the hosts. Comparative analysis with 37 Agaricomycetes genomes revealed that IstB-like domains (PF01695) were conserved across taxonomically diverse mycorrhizal genomes, where the T. matsutake genome contained four copies of this domain. Three of the IstB-like genes were overexpressed in the hyphae. Similar to other ectomycorrhizal genomes, the CAZyme gene set was reduced in T. matsutake, including losses in the glycoside hydrolase genes. The T. matsutake genome sequence provides insight into the causes and consequences of genome size inflation.