Next Article in Journal
COVID-19 Associated Pulmonary Aspergillosis (CAPA): Hospital or Home Environment as a Source of Life-Threatening Aspergillus fumigatus Infection?
Previous Article in Journal
Clinical Presentations, Cluster Analysis and Laboratory-Based Investigation of Aspergillus Otomycosis—A Single Center Experience
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JoF
Volume 8
Issue 3
10.3390/jof8030314
jof-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chemical Constituents and Molecular Mechanism of the Yellow Phenotype of Yellow Mushroom (Floccularia luteovirens)

by
Xiaolong Gan
1,2,3,†,
Xuemei Bao
1,3,4,†,
Baolong Liu
1,3,†,
Yun Li
1,3,
Dong Cao
1,3,
Hg Zhang
1,3,* and
Yuan Zong
1,2,*
1
Qinghai Province Key Laboratory of Crop Molecular Breeding, Xining 810008, China
2
Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
School of Education, Qinghai Normal University, Xining 810008, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
J. Fungi 2022, 8(3), 314; https://doi.org/10.3390/jof8030314
Submission received: 22 February 2022 / Revised: 15 March 2022 / Accepted: 17 March 2022 / Published: 18 March 2022
(This article belongs to the Section Fungal Genomics, Genetics and Molecular Biology)
Browse Figures
Review Reports Versions Notes

Abstract

:
(1) Background: Yellow mushroom (Floccularia luteovirens) is a natural resource that is highly nutritional, has a high economic value, and is found in Northwest China. Despite its value, the chemical and molecular mechanisms of yellow phenotype formation are still unclear. (2) Methods: This study uses the combined analysis of transcriptome and metabolome to explain the molecular mechanism of the formation of yellow mushroom. Subcellular localization and transgene overexpression techniques were used to verify the function of the candidate gene. (3) Results: 112 compounds had a higher expression in yellow mushroom; riboflavin was the ninth most-expressed compound. HPLC showed that a key target peak at 23.128 min under visible light at 444 nm was Vb2. All proteins exhibited the closest relationship with Agaricus bisporus var. bisporus H97. One riboflavin transporter, CL911.Contig3_All (FlMCH5), was highly expressed in yellow mushrooms with a different value (log2 fold change) of −12.98, whereas it was not detected in white mushrooms. FlMCH5 was homologous to the riboflavin transporter MCH5 or MFS transporter in other strains, and the FlMCH5-GFP fusion protein was mainly located in the cell membrane. Overexpression of FlMCH5 in tobacco increased the content of riboflavin in three transgenic plants to 26 μg/g, 26.52 μg/g, and 36.94 μg/g, respectively. (4) Conclusions: In this study, it is clear that riboflavin is the main coloring compound of yellow mushrooms, and FlMCH5 is the key transport regulatory gene that produces the yellow phenotype.

Graphical Abstract

1. Introduction

Floccularia luteovirens (yellow mushroom), also known as Armillaria luteovirens, Tricholoma luteovirens, and golden mushrooms, belongs to the genus Tricholoma of phylum Basidiomycota and family Agaricaceae, which mainly grows in Qinghai, Sichuan, Tibet and other high-altitude areas [1,2]. Yellow mushroom is expensive due to its high nutrient contents, delicious taste, and scarcity. Altitude, season, environment, and other conditions limit the domestication and cultivation of yellow mushrooms (Floccularia luteovirens) [1]. Once the growth environment of yellow mushrooms is artificially damaged, it is difficult to repair [3]. At present, research on yellow mushroom is mainly focused on its nutritional characteristics, its effects on grassland vegetation and soil, its fermentation conditions, and the determination of trace elements and fatty acids [4,5,6,7]. It has been reported that the contents of eight kinds of essential amino acids, particularly lysine and arginine, are higher in yellow mushroom compared with those in Pleurotus eryngii and Flammulina velutipes, which are frequently consumed by humans [8,9]. In addition, yellow mushroom also contains 17 minerals, with K, P, and S being predominant, followed by Fe, Zn, and Cu. Se enrichment also increases the nutritional value of yellow mushroom [2,10]. As it is a food with high nutritional value, research on yellow mushroom is imminent, but the molecular mechanisms of the yellow phenotype are still unclear. So far, a number of studies have reported that carotenoids, flavonoids, and other compounds contribute to the yellow flower color or fruit color of angiosperms, but few studies focus on fungi [11,12]. In the absence of chloroplasts, the yellow phenotype of yellow mushroom may be quite different [13,14].
Riboflavin, also known as vitamin B2, is one of the water-soluble vitamins necessary for the human body, and its deficiency can lead to oral ulcers [15,16]. The human body itself cannot synthesize riboflavin and needs to obtain it from the outside world [17]. Riboflavin exists in various organs of the human body, promoting metabolism and normal vision [18]. Thus, riboflavin plays an important role in cell proliferation and growth [19]. Riboflavin is an isopropylamine derivative with a side chain of ribitol, mainly formed from the direct precursors GTP and ribulose 5-phosphate through a series of reactions catalyzed by riboflavin synthase [20]. GTP cyclohydrolase II catalyzes the ring-opening hydrolysis reaction of GTP imidazole to generate dopamine and adenosine 3′5′-monophosphate-regulated phosphor-protein (DARPP) [21]. Riboflavin synthase β-subunit catalyzes the conversion of DHBP and 5-amino-6-ribol amino-2,4 (1H, 3H)-pyrimidine Dione (ArP) to 6,7-dimethyl-8-ribitol-2,4-dioxatetrahydrobiopterin (DRL). Finally, DRL is directly converted to riboflavin [22]. Two molecules of DRL undergo a disproportionation reaction under the action of riboflavin synthase α subunit to produce one molecule of ArPP and one molecule of riboflavin [21,23].
Among the many membrane proteins, a large protein of transporters, known as the major facilitator superfamily (MFS), has been found to have 12 transmembrane helix domains [24]. MCH5 was homologous to mammalian monocarboxylate transporters, and the MFS structural domain it contains is the main reason for its transport function. Transport experiments in S. cerevisiae and Schizosaccharomyces pombe have revealed that Mch5p is a high-affinity transporter (Km = 17 μm) with an optimal pH of 7.5 and that the expression of MCH5 is regulated by the riboflavin content in cells [25]. This result indicates that S. cerevisiae possesses a mechanism for sensing riboflavin and avoiding riboflavin deficiency by increasing the expression of the plasma membrane transporter MCH5 [26]. Additionally, MFS transporters also occupy a very important position in plants, including plant P transporters, which mediate the absorption and transport of P by roots and maintain the steady state of P in plants [27,28,29].
In the present study, metabolomic analysis was used to understand the differences between the main compounds in yellow mushroom (Floccularia luteovirens) and white mushroom (Agaricus campestris), as well as the major chemical constituents responsible for the yellow phenotype. At the same time, the major yellow compound was separated by HPLC. Second-generation transcriptome sequencing can quantify gene expression, mine transcripts, and provide reliable information to study molecular anabolic pathways in plants. Transcriptome-sequencing technology was employed to determine the major differences in the riboflavin anabolic pathway in the formation of the yellow phenotype of mushrooms and screen out the key gene that regulates the generation of the yellow and white differential phenotypes of Floccularia luteovirens.

2. Results

2.1. Analysis of Chemical Contents in Yellow Mushrooms

Overlapping display analysis of total ion chromatogram (TIC) for MS detection and analysis of different quality-control samples revealed sufficient reproducibility and reliability of the data for further investigation (Figure S1). A total of 573 compounds from 22 classes were detected (Table S1), among which 227 compounds from 20 classes exhibited significantly different levels of accumulation in the yellow and white mushrooms (Table S2). The contents of 112 compounds were higher in yellow mushroom (Figure 1a). In addition, one quinone and five steroids were detected in both yellow and white mushrooms, which did not exhibit differential accumulation (Table 1). The contents of 11 classes of compounds in yellow mushrooms were higher than those in white mushrooms. Among them, flavone, flavonol, isoflavone, carbohydrates, terpene, vitamins and their derivatives, and lipids were more than two-fold higher in yellow mushrooms, with flavone presenting more than eight-fold higher content compared with that in white mushrooms (Table 1). In addition, the contents of both flavonol and vitamins and their derivatives were significantly different between the two species of mushrooms. Three differentially accumulated flavonol compounds, dihydromyricetin and ayanin, were only found in the white mushrooms (Table S2). A relatively high content of alcohol species was observed in white mushrooms, whereas little or no alcohol content was noted in yellow mushrooms with a log2 fold change between −3 and −17. In contrast, the contents of alkaloids and indole derivatives were higher in white mushrooms than those in yellow mushrooms.
Curcumin, a differentially expressed yellow-colored polyphenol, was only found in yellow mushrooms, with the log2 fold change value reaching 11.43, although its quantitative value was low. It was found that the content of lysoPC 16:0 in yellow mushrooms was the highest, followed by galactinol and azelaic acid, but they were not yellow compounds. When sorting and screening with the peak area as the estimated content value, it was found that riboflavin held ninth place in terms of content and was dominant in both mushroom species. The riboflavin content of yellow mushroom was high, with the log2 fold change value reaching 1.53.

2.2. Identification and Content Determination of Riboflavin

After the water extraction of 5 g of fresh yellow mushrooms, HPLC results showed a key target peak at 23.128 min under visible light at 444 nm. Subsequent HPLC detection revealed that the yellow water and riboflavin standard (Vb2) were the same, with a peak at 23.126 min (Figure 1b). Quantitative detection of riboflavin in white and yellow mushrooms using the national standard method showed that the white mushrooms extract appeared dark yellow after water extraction, whereas the yellow mushrooms changed from yellow to white after water extraction, with the extract appearing bright yellow. The white and yellow mushrooms showed consistent peaks at 444 nm; retention times of 7.935 and 7.934 min; peak areas of 10,262,820 and 16,221,615; and riboflavin contents of 0.89 and 1.35 mg/100 g, respectively (Figure 1c). The results of riboflavin content were consistent with the difference in the riboflavin content in metabolome sequencing.

2.3. Transcriptome Sequencing and Data Analysis

The total number of bases obtained for each sample after transcriptome sequencing was 8.31 Gb, which reduced to 7.54 Gb after filtering (Table S3). The clean reads were further assembled into 48,997 unigenes, and after removing redundant genes, 43,750 genes were obtained. A BlastX search attained a total of 43,750 predicted proteins (Table S4). Protein databases such as Nr, Nt, Swissprot, KEGG, KOG, Interpro, GO, and Intersection were used to annotate 35,796, 13,529, 21,680, 24,808, 21,199, 24,460, 6899, and 2348 genes, respectively (Figure 2a). Homology comparison revealed that all the proteins in the mushroom samples exhibited the closest relationship with Agaricus bisporus var. bisporus H97 (20.32%), followed by Agaricus bisporus var. burnettii JB137-S8 (17.25%) (Table S5).
A total of 3863 unigenes were downregulated, and 2309 unigenes were upregulated in white mushrooms (Figure 2b). Through KEGG annotation, all the differentially expressed genes (DEGs) were annotated into 21 types of pathways and could be divided into 5 categories (Figure S5). Subsequently, 11 structural genes involved in riboflavin biosynthesis and metabolism and 60 homologous genes were selected. Among them, five structural genes were not differentially expressed, or their expression levels were low in yellow mushrooms (Table S7). The rest of the six structural genes were upregulated at different levels in yellow mushrooms, of which the transcripts of three structural genes, rib3, rib4, and FLAD1, were only found in yellow mushrooms but not in white mushrooms (Figure 2c). Furthermore, the transcript of the homologous gene of riboflavin transporter CL911.Contig3_All (FlMCH5) was highly expressed in yellow mushrooms, with a log2 fold change value of −12.98, whereas it was not detected in white mushrooms (Table S8). Therefore, FlMCH5 was screened as the key regulatory gene of the yellow phenotype for function-prediction analysis.

2.4. Conjoint Analysis of Transcriptome and Metabolome

The riboflavin cell transport pathway showed that riboflavin was upregulated by 1.54-fold in yellow mushrooms, whereas the downstream product FAD was downregulated by 1.25-fold, and the resultant glutamic acid was upregulated 1.01-fold (Figure 2c). To further evaluate the effects of transcriptome change on the metabolome, the DEGs and differentially accumulated compounds were investigated using the KEGG database. A total of 66 pathways were determined to possess both DEGs and differentially accumulated compounds, and the p value was used as the limit to determine the enrichment of genes or compounds. DEGs were mainly enriched in seven metabolic pathways, namely, Aminoacyl-tRNA Biosynthesis, Yeast Cell Cycle, Cysteine and Methionine Metabolism, Longevity-Regulating Pathway—Multiple Species, Lysine Degradation, Propanoate Metabolism, and Riboflavin Metabolism. Differentially accumulated compounds were mainly enriched in two metabolic pathways, namely, Glycerophospholipid Metabolism and Purine Metabolism (Figure S6). The CL911.Contig3_All gene was involved in Glycine, Serine, and Threonine Metabolism; Beta-Alanine Metabolism; Phenylalanine Metabolism; and Tyrosine Metabolism pathways (Table S9).

2.5. Molecular Characteristics of FlMCH5

Phylogenetic tree analysis revealed that the CL911.Contig3_All gene was homologous to the riboflavin transporter MCH5 or MFS transporter in other species (Figure 3a). Bioinformatics analysis showed that the 1–301 amino acid sequence of the CL911.Contig3_All gene belongs to the MFS, with structural domains similar to CfMFS (Crassisporium funariophilum) and HmMCH5 (Hypsizygus marmoreus) (Figure 3b). Among the CL911.Contig3_All gene protein sequences, the widest coverage was exhibited by the 7bp3.1.A protein, which belonged to the MFS and was a homodimer. These results confirmed that the CL911.Contig3_All gene was very likely to participate in importing riboflavin in yellow mushrooms.
The FlMCH5-GFP fusion cistron was constructed from cauliflower mosaic virus 35S promoter (pK7WGF-2) and was transiently expressed in tobacco cells. The results showed that the FlMCH5-GFP fusion protein was mainly located in the cell membrane, while the control GFP protein was distributed throughout the cell (Figure 3c).

2.6. Overexpression of FlMCH5 Induces Riboflavin in Tobacco

The PC2300: FlMCH5 plasmid was transferred into Agrobacterium strain GV3101 using the freeze–thaw method. The Agrobacterium-mediated leaf disk transformation method was used to obtain transgenic tobacco. For further experiments, the T3 family lines carrying the objective gene without separation were used. The color of the transgenic tobacco leaf extract is darker than the wild type. Using Vb2 standard products to calculate the content of riboflavin in tobacco, we obtained totals of 4.25 μg/g (WT), 26 μg/g (FlMCH5-1), 26.52 μg/g (FlMCH5-2), and 36.94 μg/g (FlMCH5-3) (Figure 3d). The relative riboflavin concentration of the transgenic lines was much higher than that of the wild type. These results showed that FlMCH5 could activate riboflavin biosynthesis by acting as an MFS factor in tobacco.

3. Discussion

Currently, research on yellow mushrooms is mainly focused on the determination and analysis of trace elements, whereas epigenetic investigation is still scarce. In the present study, chemical analysis suggested that the yellow phenotype of yellow mushrooms is the result of riboflavin accumulation. Stepwise MIM-EPI and RNA-Seq were confirmed to be effective strategies for performing comprehensive metabolomic and transcriptomic analyses of yellow and white mushrooms, respectively. A total of 573 compounds from 22 classes were detected, and 43,750 unigenes were assembled.

3.1. Metabolome Clarify the Main Compounds in Mushrooms

Principal Component Analysis (PCA) showed significant differences between the yellow and white mushrooms (Figure S2). OPLS-DA revealed that the data point was farthest from the origin, indicating that the greater the contribution of this point to separation between the groups, the better the correlation. The vast majority of the distances shown in the figure are red, upregulated compounds, signifying that the yellow phenotype of yellow mushrooms is mainly caused by the upregulated compounds (Figure S3). The reason for the significant difference in the flavone content between the two varieties of mushrooms was mainly the presence of sakuranetin, C-hexosyl-luteolin O-hexosyl-O-pentoside, tricin 7-O-hexosyl-O-hexoside, and tricin 5-O-hexosyl-O-hexoside in yellow mushrooms, which were almost undetectable in white mushrooms. The metabolomic analysis also found some compounds in yellow mushrooms, such as vanillic acid, 6-hydroxy-4-methylcoumarin, skimmine, curdione, α-ionone, phytocassane D, and eudesmic acid, that were not detected in white mushrooms (Table S2). These compounds have high medicinal values, including antibacterial, antitumor, anti-inflammatory, and antioxidant properties and enhance the health benefits of yellow mushrooms [30,31,32,33,34].
The highest order of compound content is mainly lipids, organic acids and derivatives, and vitamins and derivatives, but among these compounds, only riboflavin is yellow, and the other compound monomers are white powder. Therefore, the high riboflavin content is likely to be the chromogenic substance of the yellow mushroom phenotype. The quantification results of riboflavin in Cantharellus were consistent with those of the metabolomic analysis, revealing that the riboflavin content was higher in the yellow mushrooms. Interestingly, riboflavin could also be extracted from white mushrooms. The extract appeared yellow in color, further suggesting that riboflavin in white mushrooms is likely to be inside the cell and not transported outside the cell, resulting in sufficient riboflavin content but not the yellow phenotype in white mushrooms.

3.2. Riboflavin Is the Main Chemical Basis of the Yellow Phenotype

The wavelength of visible light is in the range of 400–760 nm. In nature, most of the major yellow compounds in yellow flowers, fruits, or tubers are carotenoids [35]. If an object absorbs violet light at a wavelength of 400–435 nm, then it will appear yellow-green in color [36]. The absorption ratio of chlorophyll and carotenoid is the highest at 400–520 nm (blue light), which has the maximum impact on photosynthesis [37]. Therefore, 444 nm was selected for scanning the compound that causes the yellow phenotype of yellow mushroom. A white phenotype was observed after soaking the yellow mushrooms in water, and the extract obtained was yellow in color and translucent, indicating that the compounds that cause the yellow phenotype were easily soluble in water (Figure 1c). HPLC analysis of the yellow mushroom extract at 444 nm revealed two main absorption peaks at 2.8 and 23.128 min, with the first 5 min absorption peak considered as the solvent peak. Therefore, it was speculated that the yellow phenotype might have been caused by the compound that appeared at 23.126 min. Although there were other sensitive peaks at 5–10 min intervals, the sensitive voltage value was lower. Subsequently, the compound that appeared at 23.126 min was eluted using acetonitrile, concentrated, and dried to obtain a yellow powder. Previously, 1H-NMR and 13C-NMR analysis confirmed the molecule of riboflavin extract was consistent with VB2(C17H20N4O6). This is enough to show that the yellow compound in yellow mushroom is easily soluble in water, and the only compound that appears under the yellow wavelength absorption peak is riboflavin.

3.3. Riboflavin Synthesis Also Exists in White Mushroom

Differentially accumulated compounds were mainly enriched in two metabolic pathways, namely, Glycerophospholipid Metabolism and Purine Metabolism. The KEGG enrichment map shows that Purine Metabolism has the most differential metabolic compounds (Figure S4). The Purine Metabolism downstream product is also the riboflavin synthesis substrate GTP [38,39]. In the riboflavin biosynthetic metabolic pathway, only rib2 increased by 1.88-fold; the intermediate product 5-amino-6-ribitylamino-2,4-(1H,3H)-pyrimidinedione-5′-phosphate was not differentially expressed. Furthermore, the pentose phosphate pathway was finally regulated by rib3, rib4, and ribE to form riboflavin. Interestingly, rib3, rib4, RFK, and FLAD1 were upregulated in yellow mushrooms. RibE was only expressed in white mushrooms but not in yellow mushrooms. The specific expression of ribE in white mushrooms might be the main reason for the extraction of riboflavin from white mushrooms, although the content of extracted riboflavin was not as high as that from yellow mushrooms [40,41]. The regulation of structural genes related to riboflavin metabolism can only explain the low riboflavin content in white mushrooms and cannot clarify the yellow phenotype; hence, it is likely that the yellow phenotype might be closely related to cell transport.

3.4. FlMCH5 Produces the Yellow Phenotype

In the early stages, the genome assembly of yellow mushroom was completed, and the genome of yellow mushroom is 28.7 Mb, with BUSCO reaching 93.9% [42]. Evolutionary analysis showed that yellow mushroom was near Agaricus bisporus, and FlMCH5 could also be retrieved in the yellow mushroom genome library, but it did not exist in Agaricus bisporus or white mushroom in this study.
In S. cerevisiae, proline degradation has been reported to be closely associated with riboflavin transport [43,44], while only glycyl-L-proline was downregulated in the metabolome of yellow mushrooms, with a log2 fold change value reaching −1.32. Furthermore, in S. cerevisiae, proline degradation can activate the PUT3 transcription factor [45]. PUT3 regulates the other two proline-regulating genes, PUT1 and PUT2, and can also positively control the riboflavin transporter MCH5 [46]. In the present study, transcriptome comparison revealed that the homologous genes of PUT3 and FlMCH5 in yellow mushrooms were Unigene2188_All and CL911.Contig3_All, respectively, and that these two transcripts were only expressed in yellow mushrooms and not in white mushrooms. While PUT2 homologous genes could be detected in all transcripts, PUT1 homologous genes could not be determined, and the expression increase in PUT2 in white mushrooms was only 0.14. Western blot analyses showed that the increased riboflavin expression in the mutant strain was due to the activation of MCH5 mediated by PUT3. The gel shift test proved that PUT3 can bind to the promoter of MCH5. It must be noted that PUT3-mediated transcriptional activation requires proline as an inducer [47]. Thus, downregulation of proline in chanterelles resulted in the upregulation of the transcription factor PUT3 (Unigene2188_All), which activated the riboflavin transporter FlMCH5 (CL911.Contig3_All) that transported extracellular riboflavin into the cell, leading to the yellow phenotype (Table S7).

3.5. FlMCH5 Is a Functional MFS Transcription Factor Gene Regulation the Yellow Phenotype

FlMCH5 belongs to MFS and was homologous to the riboflavin transporter MCH5 in Hypsizygus marmoreus. MFS is one of the largest families of transporters. Analysis of the amino acid structure of 66 genes revealed that the gene with the widest coverage was 7bp3.1.A, which had a structure that conformed to the condition that each domain of the MFS protein was composed of inverted 3 + 3 repeats [47]. Subcellular localization showed that FlMCH5-GFP acted more on the cell wall, further proving that the gene was related to membrane transport. Although there is no phenotypic change in transgenic tobacco, there are significant differences in riboflavin extract in leaves. The color of the leaf extract of transgenic tobacco is deeper, and the riboflavin content is higher. All of these results implied that FlMCH5 was a functional MFS transcription factor regulating riboflavin biosynthesis.

3.6. Creation of a New Disease-Resistant Tobacco with a High Riboflavin Content

Riboflavin in tobacco was discovered as early as the late 1950s; as the elicitor of tobacco defense response, riboflavin can enhance the resistance of tobacco to diseases such as fungi, bacteria, and viruses. At the same time, riboflavin in tobacco can be used as an antioxidant, which helps prevent the conversion of nitrate in the smoke, and can also be used as a photosensitizer to catalyze the oxidation and degradation of nicotine by antibodies. The content of riboflavin in tobacco from Guizhou, Sichuan, and other places in China determined by LC-ITMS was 5.81, 5.77, and 8.50 μg/g, respectively. In this study, the riboflavin content in tobacco overexpressing the FlMCH5 gene was as high as 26 μg/g–36.94 μg/g. Compared with wild-type or other domestic tobacco, the riboflavin content increased by 5–7 times. Subsequent cross-breeding work has a high probability of cultivating a new type of tobacco with strong disease resistance.

4. Materials and Methods

4.1. Plant Materials

The yellow mushroom (Floccularia luteovirens) (nwipbYM-1807) and white mushroom (Agaricus campestris) (nwipbYM-1808) were collected from the Qilian Mountains in Qinghai Province, China (latitude: 39°24′ N, longitude: 98°53′ E, altitude approximately 3800 m). Two kinds of mushrooms were identified by morphological identification as Floccularia luteovirens and Agaricus campestris, and the collected samples were stored at Northwest Plateau Institute of Biology, Chinese Academy of Sciences, China. White mushroom is a white mutant within 100 m of the colony circle of yellow mushroom. No permission was required for sample collection.

4.2. Analysis of Chemical Contents

In this study, white mushroom and yellow mushroom were selected for genetic analysis under the same environmental conditions. The white and yellow mushroom samples were frozen in liquid nitrogen immediately after collection. Before chemical and transcriptomic analysis, the white and yellow caps were stored in a refrigerator at −80 °C. The lyophilized sample was crushed at 30 Hz for 1 min using a steel ball crushing mill (Retsch MM 400, Haan, Germany). Then, 100 mg of the resultant powder was extracted with 1.0 mL of pure water overnight at 4 °C in dark. Then, the solution was centrifuged at 12,000× g for 10 min, and the extract was absorbed (CNWBOND Carbon-GCB SPE column, 250 mg, 3 mL; ANPEL, Shanghai, China, www.anpel.com.cn/cnw (accessed on 28 June 2019)), filtered (SCAA-104, pore size, 0.22 μm; ANPEL, Shanghai, China, http://www.anpel.com.cn/ (accessed on 28 June 2019), and subjected to liquid chromatography–mass spectrometry (LC/MS) analysis. Then, the extracts were analyzed using LC electrospray tandem MS (ESI MS/MS) system (Applied Biosystems 4500 QTRAP, http://www.appliedbiosystems.com.cn/ (accessed on 29 June 2019)) with water (0.01% acetic acid) and acetonitrile (0.01% acetic acid) as mobile phases [48]. The gradient program was as follows: 5% acetonitrile at 0 min, 20% acetonitrile at 3 min, 30% acetonitrile at 8 min, and 90% acetonitrile at 13 min; flow rate of 0.40 mL/min; temperature of 40 °C; and injection volume of 2 μL.
The raw data of LC/MS were converted into CDF format (NetCDF) using Agilent 5975c data analysis software and then preprocessed by XCMS (www.bioconductor.org (accessed on 8 August 2019)) to export the three-dimensional data matrix in CSV format [49,50]. Unsupervised principal components (PCA) were used to study the overall distribution between samples and stability of the entire analysis process. Subsequently, supervised partial least squares analysis (PLS-DA) was employed to distinguish the metabolic profiles between the groups. The model was effective if the overall difference and R2 intercept were less than 0.4 and Q2 intercept was less than 0.05; otherwise, it was defined as simulation overfitting. After determining the model, the difference in metabolites between the groups was analyzed. In the PLS-DA, the variable weight value (VIP) was higher than 1, and the variable that was less than 0.05 through t-test was considered to be the difference variable. To prevent the PLS-DA model from overfitting, the method of seven cycles of interactive verification and 999 tests of response ranking was used to examine the quality of the model [51,52]. R (www.r-project.org/ (accessed on 23 August 2019)) thermogram cluster analysis was used to obtain more obvious metabolite distribution [53].

4.3. Qualitative and Quantitative Analysis of Yellow Compound

Using a vacuum freeze dryer, the lyophilized sample was weighed, ground into a crushed residue, mixed with adequate amount of sterile water, and sonicated at 50 °C for 30 min in a water bath. Then, the supernatant yellow extract was pipetted, filtered through a 0.45 μm membrane, and subjected to HPLC under the following conditions: 0–5 min, 5–20% acetonitrile; 5–30 min, 20–25% acetonitrile; 30–35 min, 25–90% acetonitrile; column temperature, 25 °C; injection volume, 50 μL; operation time, 35 min; and spectral scanning value, 444 nm [54]. The target peak was selected for preparation and recovery. The injection volume was 1.8 mL, number of injections was 50 times, and detection wavelength was 444 nm. The prepared solution was concentrated and freeze-dried to obtain a dry powder of riboflavin extract [55]. A total of 0.05 g of riboflavin standard (Vb2: Yuanye Biology, CAS#83-88-5, HPLC ≥ 99%, experimental concentration 1%) was dissolved in 2 mL methanol for preparation liquid check. The riboflavin content in the yellow and white mushrooms was determined according to national standards (GB5009.85-2016). The concentration of Vb2 standard stock solution was 2.5 mg/g; we diluted the concentration to 0.02, 0.04, 0.06, and 0.08 mg/g; and the injection volume was 10 µL. Y = 4570.3X and R2 = 0.9989.

4.4. Transcriptome Sequencing and Analysis

The white and yellow mushroom caps were frozen at −80 °C, and the total RNA was extracted according to sequencing requirements. Mushrooms RNA was extracted by using TIANGEN polysaccharide and polyphenol extraction kit (CAT#:DP441, Beijing, China. TIANGEN Biotech Company). A total of 1 μL of extracted RNA solution was aspirated, and the RNA concentration was determined using NanoDrop. The ratio of A260/A280 was about 1.8–2.1. The integrity and high purity of the extracted RNA were determined using NanoDrop, Qubit2.0, and Agilent 2100 Bioanalyzer [56]. After qualitative analysis of the samples, the sequencing library was constructed. The selected cDNA was sequenced using Illumina high-throughput sequencer (X-10) [57]. After obtaining high-quality sequencing data, Trinity was employed for sequence assembly. The k-mer library was constructed by interrupting the reads according to the short k-mer fragment (k-mer) [58]. We used HISAT2 v2.1.0 to align the quality-controlled transcriptome reads to the yellow mushroom (Floccularia luteovirens) genome (https://doi.org/10.6084/m9.figshare.19354028 (accessed on 15 October 2019)).
BLAST software was utilized to compare the Unigene sequence with NR, Swiss-Prot, GO, COG, KOG, and KEGG databases; HMMER software was employed to compare the Unigene sequence with Pfam database after predicting the amino acid sequence of Unigene; and the annotation information of Unigene was obtained [59]. By selecting BLAST parameter E-value less than 1 × 10−5 and HMMER parameter E-value less than 1 × 10−10, Unigene with annotation information was finally acquired. Bowtie was employed to compare the sequenced reads with the Unigene library and combined with RSEM to evaluate the expression level. FPKM value was used to indicate the expression abundance of the corresponding Unigene. As the identification standard for the screening process, the FDR was less than 0.001, and the difference in expression between the two groups’ fold change was higher than 4 [60].

4.5. Functional Analysis of FlMCH5

The amino acid structure was analyzed using transmembrane transporter (https://swissmodel.expasy.org/interactive (accessed on 7 March 2020)), and the protein was identified (https://www.rcsb.org/structure (accessed on 7 March 2020)). The MEGA7 software was used for the development of tree drawings [61]. The “TAG” sequence of the FlMCH5 gene was removed and constructed in a pK7WGF-2 expression vector with the GUS gene by using gateway technology. After shaking Agrobacterium overnight, the supernatant was centrifuged to remove the supernatant, and the permeate was added. After diluting 10 times, the OD600 value was measured and controlled at about 0.4, and it was allowed to stand at room temperature for 1–3 h. We put the tobacco under a white fluorescent lamp for 1 h, injected the suspended bacteria liquid into the back of the tobacco using a syringe, and cultivated it in the dark for 2–3 days after the injection. We cut off the infected area and observed with green fluorescence after filming [62].

4.6. Overexpression of FlMCH5 in Tobacco

The full length of FlMCH5 was cloned by PCR; the primers for the selected genes were designed by Vector NTI10.0 (Table S6) [63]. The sequence was cloned into the pEASY®-Blunt vector (TransGen Biotech, Beijing, China), which transformed into Escherichia coli. DH5α cell and the positive cloned were sequenced. FlMCH5 plasmid and PC2300s plasmid were double-digested using XbaI and BamHI. The PC2300s:FlMCH5 plasmid was transformed into GV3101. Nicotiana tabacum (Samsun) was used for transformation using the leaf disc transformation method [64]. Quantitative analysis of riboflavin in transgenic plants was performed using the same method as for yellow mushrooms.

5. Conclusions

This study used the combined analysis of the transcriptome and metabolome to explain the molecular mechanism of the formation of yellow mushroom, a characteristic fungus resource in Northwest China. We found that the formation of the yellow phenotype was strongly associated the accumulation of riboflavin and the regulation of the transporter FlMCH5, which transports riboflavin from outside the cell into the cell, thus creating a yellow phenotype.
This study not only confirmed that yellow mushroom has a variety of amino acids but also found that Vb2(riboflavin) laid a sloid foundation for the yellow phenotype. Our research provides a new opportunity for advertising the yellow mushroom in Northwest China. In addition to yeast, tobacco can also be used as one of the transporter function analysis carriers.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof8030314/s1, Figure S1. TIC overlaps by QC sample mass spectrometry. Figure S2. PCA cluster analysis of the metabolome. Figure S3. OPLS-DS cluster analysis of compounds. Figure S4. All genes’ KEGG pathway enrichment. Figure S5. KEGG pathway enrichment of differentially expressed genes. Figure S6. Combined enrichment analysis of all differentially expressed genes and compounds. Table S1. All chemical compounds in yellow and white mushrooms. Table S2. The different accumulation of chemical compounds in yellow and white mushrooms. Table S3. Distribution of the clean reads of the sequencing. Table S4. Distribution of the lengths of the assembled unigenes. Table S5. Phylogenetic analysis of KOG species. Table S6. The primers used in this study. Table S7. The expression differences of structural and regulatory genes in the riboflavin biosynthesis pathway. Table S8. The different expressions of unigenes in yellow and white mushrooms. Table S9. The KEGG pathways of different expressions of unigenes and different accumulations of chemical compounds.

Author Contributions

Conceptualization, Y.Z. and H.Z.; methodology, B.L.; software, X.G.; validation, X.B. and B.L.; formal analysis, B.L.; investigation, Y.L.; resources, D.C.; data curation, X.G.; writing—original draft preparation, Y.Z.; writing—review and editing, Y.Z.; visualization, B.L.; supervision, Y.L.; project administration, B.L.; funding acquisition, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Qinghai Science and Technology Department for Yuan Zong project: 2021-ZJ-955Q.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The transcriptomic data were successfully uploaded to NCBI, PRJNA647230: RNA-seq data of Agaricus campestris, PRJNA647227: RNA-seq data of Floccularia luteovirens. All data generated or analyzed during this study are included within the article and its additional files.

Acknowledgments

We are grateful for a project funded by the “Western Light Talent Program” of the Chinese Academy of Science.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xing, R.; Gao, Q.B.; Zhang, F.Q.; Fu, P.C.; Wang, J.L.; Yan, H.Y.; Chen, S.L. Genetic variation and phylogenetic relationships of the ectomycorrhizal Floccularia luteovirens on the Qinghai-Tibet Plateau. J. Microbiol. 2017, 55, 600–606. [Google Scholar] [CrossRef]
  2. Liu, Z.; Lu, H.; Zhang, X.; Chen, Q. 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. J. Fungi 2021, 7, 887. [Google Scholar] [CrossRef]
  3. Liu, Z.; Jiao, Y.; Lu, H.; Shu, X.; Chen, Q. Chemical characterization, antioxidant properties and anticancer activity of exopolysaccharides from Floccularia luteovirens. Carbohydr. Polym. 2020, 229, 115432. [Google Scholar] [CrossRef]
  4. Fu, Y.; Yu, Y.; Tan, H.; Wang, B.; Peng, W.; Sun, Q. Metabolomics reveals dopa melanin involved in the enzymatic browning of the yellow cultivars of East Asian golden needle mushroom (Flammulina filiformis). Food Chem. 2022, 370, 131295. [Google Scholar] [CrossRef]
  5. Ng, T.B.; Wang, H.X. A homodimeric laccase with unique characteristics from the yellow mushroom Cantharellus cibarius. Biochem Biophys Res. Commun. 2004, 313, 37–41. [Google Scholar] [CrossRef]
  6. Vinogradov, E.; Petersen, B.O.; Duus, J.O.; Wasser, S. The structure of the glucuronoxylomannan produced by culinary-medicinal yellow brain mushroom (Tremella mesenterica Ritz.:Fr., Heterobasidiomycetes) grown as one cell biomass in submerged culture. Carbohydr. Res. 2004, 339, 1483–1489. [Google Scholar] [CrossRef]
  7. Wang, H.X.; Ngai, H.K.; Ng, T.B. A ubiquitin-like peptide with ribonuclease activity against various polyhomoribonucleotides from the yellow mushroom Cantharellus cibarius. Peptides 2003, 24, 509–513. [Google Scholar] [CrossRef]
  8. Xing, R.; Yan, H.Y.; Gao, Q.B.; Zhang, F.Q.; Wang, J.L.; Chen, S.L. Microbial communities inhabiting the fairy ring of Floccularia luteovirens and isolation of potential mycorrhiza helper bacteria. J. Basic Microbiol. 2018, 58, 554–563. [Google Scholar] [CrossRef] [PubMed]
  9. Rosnina, A.G.; Tan, Y.S.; Abdullah, N.; Vikineswary, S. Morphological and molecular characterization of yellow oyster mushroom, Pleurotus citrinopileatus, hybrids obtained by interspecies mating. World J. Microbiol. Biotechnol. 2016, 32, 18. [Google Scholar] [CrossRef] [PubMed]
  10. Li, Q.; Liu, J.; Shang, X.; Li, Y.; Zhang, L.; Li, Z.; Jiang, N.; Tan, Q.; Yu, H.; Song, C. Characterizing Diversity Based on the Chemical and Nutritional Composition of Shiitake Culinary-Medicinal Mushroom Lentinula edodes (Agaricomycetes) Commonly Cultivated in China. Int. J. Med. Mushrooms 2021, 23, 51–64. [Google Scholar] [CrossRef] [PubMed]
  11. Berman, J.; Sheng, Y.; Gomez Gomez, L.; Veiga, T.; Ni, X.; Farre, G.; Capell, T.; Guitian, J.; Guitian, P.; Sandmann, G.; et al. Red Anthocyanins and Yellow Carotenoids Form the Color of Orange-Flower Gentian (Gentiana lutea L. var. aurantiaca). PLoS ONE 2016, 11, e0162410. [Google Scholar] [CrossRef] [PubMed]
  12. Hoshino, A.; Mizuno, T.; Shimizu, K.; Mori, S.; Fukada-Tanaka, S.; Furukawa, K.; Ishiguro, K.; Tanaka, Y.; Iida, S. Generation of Yellow Flowers of the Japanese Morning Glory by Engineering Its Flavonoid Biosynthetic Pathway toward Aurones. Plant. Cell Physiol 2019, 60, 1871–1879. [Google Scholar] [CrossRef]
  13. Xu, Q.; Tang, C.; Wang, X.; Sun, S.; Zhao, J.; Kang, Z.; Wang, X. An effector protein of the wheat stripe rust fungus targets chloroplasts and suppresses chloroplast function. Nat. Commun. 2019, 10, 5571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Miura, E.; Kato, Y.; Matsushima, R.; Albrecht, V.; Laalami, S.; Sakamoto, W. The balance between protein synthesis and degradation in chloroplasts determines leaf variegation in Arabidopsis yellow variegated mutants. Plant. Cell 2007, 19, 1313–1328. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Suwannasom, N.; Kao, I.; Pruss, A.; Georgieva, R.; Baumler, H. Riboflavin: The Health Benefits of a Forgotten Natural Vitamin. Int. J. Mol. Sci. 2020, 21, 950. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Saedisomeolia, A.; Ashoori, M. Riboflavin in Human Health: A Review of Current Evidences. Adv. Food Nutr. Res. 2018, 83, 57–81. [Google Scholar] [CrossRef] [PubMed]
  17. Thakur, K.; Tomar, S.K.; Singh, A.K.; Mandal, S.; Arora, S. Riboflavin and health: A review of recent human research. Crit. Rev. Food Sci. Nutr. 2017, 57, 3650–3660. [Google Scholar] [CrossRef]
  18. Mosegaard, S.; Dipace, G.; Bross, P.; Carlsen, J.; Gregersen, N.; Olsen, R.K.J. Riboflavin Deficiency-Implications for General Human Health and Inborn Errors of Metabolism. Int. J. Mol. Sci. 2020, 21, 3847. [Google Scholar] [CrossRef]
  19. Lambertz, C.; Leopold, J.; Damme, K.; Vogt-Kaute, W.; Ammer, S.; Leiber, F. Effects of a riboflavin source suitable for use in organic broiler diets on performance traits and health indicators. Animal 2020, 14, 716–724. [Google Scholar] [CrossRef]
  20. Lehmann, M.; Degen, S.; Hohmann, H.P.; Wyss, M.; Bacher, A.; Schramek, N. Biosynthesis of riboflavin. Screening for an improved GTP cyclohydrolase II mutant. FEBS J. 2009, 276, 4119–4129. [Google Scholar] [CrossRef]
  21. Fassbinder, F.; Kist, M.; Bereswill, S. Structural and functional analysis of the riboflavin synthesis genes encoding GTP cyclohydrolase II (ribA), DHBP synthase (ribBA), riboflavin synthase (ribC), and riboflavin deaminase/reductase (ribD) from Helicobacter pylori strain P1. FEMS Microbiol. Lett. 2000, 191, 191–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Liauta-Teglivets, O.; Hasslacher, M.; Boretskii Iu, R.; Kohlwein, S.D.; Shavlovskii, G.M. Molecular cloning of the GTP-cyclohydrolase structural gene RIB1 of Pichia guilliermondii involved in riboflavin biosynthesis. Yeast 1995, 11, 945–952. [Google Scholar] [CrossRef] [PubMed]
  23. Herz, S.; Eberhardt, S.; Bacher, A. Biosynthesis of riboflavin in plants. The ribA gene of Arabidopsis thaliana specifies a bifunctional GTP cyclohydrolase II/3,4-dihydroxy-2-butanone 4-phosphate synthase. Phytochemistry 2000, 53, 723–731. [Google Scholar] [CrossRef]
  24. Reddy, V.S.; Shlykov, M.A.; Castillo, R.; Sun, E.I.; Saier, M.H., Jr. The major facilitator superfamily (MFS) revisited. FEBS J. 2012, 279, 2022–2035. [Google Scholar] [CrossRef] [PubMed]
  25. Reihl, P.; Stolz, J. The monocarboxylate transporter homolog Mch5p catalyzes riboflavin (vitamin B2) uptake in Saccharomyces cerevisiae. J. Biol. Chem. 2005, 280, 39809–39817. [Google Scholar] [CrossRef] [Green Version]
  26. Dulyayangkul, P.; Satapoomin, N.; Avison, M.B.; Charoenlap, N.; Vattanaviboon, P.; Mongkolsuk, S. Over-Expression of Hypochlorite Inducible Major Facilitator Superfamily (MFS) Pumps Reduces Antimicrobial Drug Susceptibility by Increasing the Production of MexXY Mediated by ArmZ in Pseudomonas aeruginosa. Front. Microbiol. 2020, 11, 592153. [Google Scholar] [CrossRef]
  27. Pasqua, M.; Grossi, M.; Zennaro, A.; Fanelli, G.; Micheli, G.; Barras, F.; Colonna, B.; Prosseda, G. The Varied Role of Efflux Pumps of the MFS Family in the Interplay of Bacteria with Animal and Plant Cells. Microorganisms 2019, 7, 285. [Google Scholar] [CrossRef] [Green Version]
  28. Vishwakarma, P.; Banerjee, A.; Pasrija, R.; Prasad, R.; Lynn, A.M. Phylogenetic and conservation analyses of MFS transporters. 3 Biotech 2018, 8, 462. [Google Scholar] [CrossRef]
  29. Liu, J.; Fu, S.; Yang, L.; Luan, M.; Zhao, F.; Luan, S.; Lan, W. Vacuolar SPX-MFS transporters are essential for phosphate adaptation in plants. Plant Signal. Behav. 2016, 11, e1213474. [Google Scholar] [CrossRef] [Green Version]
  30. Obydah, W.O.; Shaker, G.A.; Samir, S.M.; El Bassiony, S.F.; Abd El Moneim, H.A. Effect of vanillic acid and exercise training on fatty liver and insulin resistance in rats: Possible role of fibroblast growth factor 21 and autophagy. Physiol. Int. 2021, 108, 412–426. [Google Scholar] [CrossRef]
  31. Popova, S.A.; Shevchenko, O.G.; Chukicheva, I.Y. Synthesis of new coumarin[1,3]oxazine derivatives of 7-hydroxy-6-isobornyl-4-methylcoumarin and their antioxidant activity. Chem. Biol. Drug Des. 2021, 1–11. [Google Scholar] [CrossRef] [PubMed]
  32. Liang, X.X.; Li, Q.; Li, H.J.; Ning, Y.M.; Zhang, R.H.; Zhang, X.J.; Li, X.L.; Xiao, W.L. Centrantheroside F, a new ionone glycoside from Centranthera grandiflora. J. Asian Nat. Prod. Res. 2021, 1–7. [Google Scholar] [CrossRef] [PubMed]
  33. Umemura, K.; Ogawa, N.; Shimura, M.; Koga, J.; Usami, H.; Kono, T. Possible role of phytocassane, rice phytoalexin, in disease resistance of rice against the blast fungus Magnaporthe grisea. Biosci. Biotechnol. Biochem. 2003, 67, 899–902. [Google Scholar] [CrossRef] [Green Version]
  34. El-Shaheny, R.; Belal, F. Green conventional and first-order derivative fluorimetry methods for determination of trimebutine and its degradation product (eudesmic acid). Emphasis on the solvent and pH effects on their emission spectral properties. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2020, 226, 117603. [Google Scholar] [CrossRef] [PubMed]
  35. Vaddamanu, M.; Sathyanarayana, A.; Masaya, Y.; Sugiyama, S.; Kazuhisa, O.; Velappan, K.; Nandeshwar, M.; Hisano, K.; Tsutsumi, O.; Prabusankar, G. Acridine N-Heterocyclic Carbene Gold(I) Compounds: Tuning from Yellow to Blue Luminescence. Chem. Asian J. 2021, 16, 521–529. [Google Scholar] [CrossRef] [PubMed]
  36. Wu, J.; Wang, X.; He, Y.; Li, J.; Ma, K.; Zhang, Y.; Li, H.; Yin, C.; Zhang, Y. Stability evaluation of gardenia yellow pigment in presence of different phenolic compounds. Food Chem. 2022, 373, 131441. [Google Scholar] [CrossRef]
  37. Mallmann, L.P.; Tischer, B.; Vizzotto, M.; Rodrigues, E.; Manfroi, V. Comprehensive identification and quantification of unexploited phenolic compounds from red and yellow araca (Psidium cattleianum Sabine) by LC-DAD-ESI-MS/MS. Food Res. Int. 2020, 131, 108978. [Google Scholar] [CrossRef]
  38. Liesegang, A.; Fritsche, W.; Schlee, D.; Reinbothe, H. [Purine metabolism and riboflavin biosynthesis in microorganisms. IV. Storing of riboflavin in the cells of Candida guilliermondii (Cast.) Lang. et G]. Zeitschrift für allgemeine Mikrobiologie. 1971, 11, 121–130. [Google Scholar] [CrossRef]
  39. Brown, E.G.; Goodwin, T.W.; Jones, O.T. Studies on the biosynthesis of riboflavin. 4. Purine metabolism and riboflavin synthesis in Eremothecium ashbyii. Biochem. J. 1958, 68, 40–49. [Google Scholar] [CrossRef] [Green Version]
  40. Fedorovych, D.V.; Dmytruk, K.V.; Sibirny, A.A. Recent Advances in Construction of the Efficient Producers of Riboflavin and Flavin Nucleotides (FMN, FAD) in the Yeast Candida famata. Methods Mol. Biol. 2021, 2280, 15–30. [Google Scholar] [CrossRef]
  41. Jin, C.; Yao, Y.; Yonezawa, A.; Imai, S.; Yoshimatsu, H.; Otani, Y.; Omura, T.; Nakagawa, S.; Nakagawa, T.; Matsubara, K. Riboflavin Transporters RFVT/SLC52A Mediate Translocation of Riboflavin, Rather than FMN or FAD, across Plasma Membrane. Biol. Pharm. Bull. 2017, 40, 1990–1995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Gan, X.; Cao, D.; Zhang, Z.; Cheng, S.; Wei, L.; Li, S.; Liu, B. Draft Genome Assembly of Floccularia luteovirens, an Edible and Symbiotic Mushroom on Qinghai-Tibet Plateau. G3 2020, 10, 1167–1173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Mukai, Y.; Kamei, Y.; Liu, X.; Jiang, S.; Sugimoto, Y.; Mat Nanyan, N.S.B.; Watanabe, D.; Takagi, H. Proline metabolism regulates replicative lifespan in the yeast Saccharomyces cerevisiae. Microb. Cell 2019, 6, 482–490. [Google Scholar] [CrossRef] [PubMed]
  44. Nugroho, R.H.; Yoshikawa, K.; Matsuda, F.; Shimizu, H. Positive effects of proline addition on the central metabolism of wild-type and lactic acid-producing Saccharomyces cerevisiae strains. Bioprocess. Biosyst. Eng. 2016, 39, 1711–1716. [Google Scholar] [CrossRef]
  45. Tebung, W.A.; Omran, R.P.; Fulton, D.L.; Morschhauser, J.; Whiteway, M. Put3 Positively Regulates Proline Utilization in Candida albicans. mSphere 2017, 2, 317–354. [Google Scholar] [CrossRef] [Green Version]
  46. des Etages, S.A.; Falvey, D.A.; Reece, R.J.; Brandriss, M.C. Functional analysis of the PUT3 transcriptional activator of the proline utilization pathway in Saccharomyces cerevisiae. Genetics 1996, 142, 1069–1082. [Google Scholar] [CrossRef]
  47. Spitzner, A.; Perzlmaier, A.F.; Geillinger, K.E.; Reihl, P.; Stolz, J. The proline-dependent transcription factor Put3 regulates the expression of the riboflavin transporter MCH5 in Saccharomyces cerevisiae. Genetics 2008, 180, 2007–2017. [Google Scholar] [CrossRef] [Green Version]
  48. Sharifi-Rad, J.; Song, S.; Ali, A.; Subbiah, V.; Taheri, Y.; Suleria, H.A.R. LC-ESI-QTOF-MS/MS characterization of phenolic compounds from Pyracantha coccinea M.Roem. and their antioxidant capacity. Cell Mol. Biol. 2021, 67, 201–211. [Google Scholar] [CrossRef]
  49. Puffeld, M.; Seiler, C.; Kuhlmann, M.; Sreenivasulu, N.; Butardo, V.M., Jr. Analysis of Developing Rice Grain Transcriptome Using the Agilent Microarray Platform. Methods Mol. Biol. 2019, 1892, 277–300. [Google Scholar] [CrossRef]
  50. Domingo-Almenara, X.; Siuzdak, G. Metabolomics Data Processing Using XCMS. Methods Mol. Biol. 2020, 2104, 11–24. [Google Scholar] [CrossRef]
  51. Saccenti, E.; Timmerman, M.E. Approaches to Sample Size Determination for Multivariate Data: Applications to PCA and PLS-DA of Omics Data. J. Proteome Res. 2016, 15, 2379–2393. [Google Scholar] [CrossRef]
  52. Worley, B.; Halouska, S.; Powers, R. Utilities for quantifying separation in PCA/PLS-DA scores plots. Anal. Biochem. 2013, 433, 102–104. [Google Scholar] [CrossRef] [Green Version]
  53. Khandakar, A.; Chowdhury, M.E.H.; Ibne Reaz, M.B.; Md Ali, S.H.; Hasan, M.A.; Kiranyaz, S.; Rahman, T.; Alfkey, R.; Bakar, A.A.A.; Malik, R.A. A machine learning model for early detection of diabetic foot using thermogram images. Comput. Biol. Med. 2021, 137, 104838. [Google Scholar] [CrossRef] [PubMed]
  54. Schwechheimer, S.K.; Becker, J.; Peyriga, L.; Portais, J.C.; Sauer, D.; Muller, R.; Hoff, B.; Haefner, S.; Schroder, H.; Zelder, O.; et al. Improved riboflavin production with Ashbya gossypii from vegetable oil based on (13)C metabolic network analysis with combined labeling analysis by GC/MS, LC/MS, 1D, and 2D NMR. Metab. Eng. 2018, 47, 357–373. [Google Scholar] [CrossRef] [PubMed]
  55. Deng, K.; Ouyang, J.; Hu, N.; Meng, J.; Su, C.; Wang, J.; Wang, H. Improved colorimetric analysis for subtle changes of powdered anthocyanins extracted from Lycium ruthenicum Murr. Food Chem. 2022, 371, 131080. [Google Scholar] [CrossRef] [PubMed]
  56. Villa-Rodriguez, E.; Ibarra-Gamez, C.; de Los Santos-Villalobos, S. Extraction of high-quality RNA from Bacillus subtilis with a lysozyme pre-treatment followed by the Trizol method. J. Microbiol. Methods 2018, 147, 14–16. [Google Scholar] [CrossRef]
  57. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 2013, 8, 1494–1512. [Google Scholar] [CrossRef]
  58. Kim, C.S.; Winn, M.D.; Sachdeva, V.; Jordan, K.E. K-mer clustering algorithm using a MapReduce framework: Application to the parallelization of the Inchworm module of Trinity. BMC Bioinform. 2017, 18, 467. [Google Scholar] [CrossRef] [Green Version]
  59. Yu, D.S.; Lee, D.H.; Kim, S.K.; Lee, C.H.; Song, J.Y.; Kong, E.B.; Kim, J.F. Algorithm for predicting functionally equivalent proteins from BLAST and HMMER searches. J. Microbiol. Biotechnol. 2012, 22, 1054–1058. [Google Scholar] [CrossRef]
  60. Zhu, D.; Hero, A.O.; Qin, Z.S.; Swaroop, A. High throughput screening of co-expressed gene pairs with controlled false discovery rate (FDR) and minimum acceptable strength (MAS). J. Comput. Biol. 2005, 12, 1029–1045. [Google Scholar] [CrossRef] [Green Version]
  61. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Zong, Y.; Xi, X.; Li, S.; Chen, W.; Zhang, B.; Liu, D.; Liu, B.; Wang, D.; Zhang, H. Allelic Variation and Transcriptional Isoforms of Wheat TaMYC1 Gene Regulating Anthocyanin Synthesis in Pericarp. Front. Plant. Sci 2017, 8, 1645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Lu, G.; Moriyama, E.N. Vector NTI, a balanced all-in-one sequence analysis suite. Brief. Bioinform. 2004, 5, 378–388. [Google Scholar] [CrossRef]
  64. Zong, Y.; Li, S.; Xi, X.; Cao, D.; Wang, Z.; Wang, R.; Liu, B. Comprehensive Influences of Overexpression of a MYB Transcriptor Regulating Anthocyanin Biosynthesis on Transcriptome and Metabolome of Tobacco Leaves. Int. J. Mol. Sci. 2019, 20, 5123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Jof 08 00314 g001 550
Figure 1. Compound analysis of yellow mushroom. (a) Scatter plot of differential accumulation compounds, where red dots represent compounds with high content in yellow mushroom. (b) Full scan chromatogram of yellow mushroom at 444 nm after ultrasonic extraction with water. The compound was prepared and freeze-dried for about 23.128 min. (c) The white and yellow mushrooms were ultrasonically extracted with water until presentation of the white phenotype, and then the riboflavin content was determined, respectively.
Figure 1. Compound analysis of yellow mushroom. (a) Scatter plot of differential accumulation compounds, where red dots represent compounds with high content in yellow mushroom. (b) Full scan chromatogram of yellow mushroom at 444 nm after ultrasonic extraction with water. The compound was prepared and freeze-dried for about 23.128 min. (c) The white and yellow mushrooms were ultrasonically extracted with water until presentation of the white phenotype, and then the riboflavin content was determined, respectively.
Jof 08 00314 g001
Jof 08 00314 g002 550
Figure 2. Comparative analysis of transcriptome between yellow and white mushrooms. (a) BLAST software was utilized to compare the Unigene sequence with NR, Swiss-Prot, GO, COG, KOG, and KEGG databases. (b) Differentially expressed genes between yellow and white mushrooms, where red represents the number of genes upregulated in white mushroom. (c) Simulation diagram of riboflavin metabolic pathway and riboflavin transport in cells. Arrows represent the direction of the composite, “+” represents that PUT3 positively regulates FlMCH5, and “–” represents that PUT3 negatively regulates PUT2.
Figure 2. Comparative analysis of transcriptome between yellow and white mushrooms. (a) BLAST software was utilized to compare the Unigene sequence with NR, Swiss-Prot, GO, COG, KOG, and KEGG databases. (b) Differentially expressed genes between yellow and white mushrooms, where red represents the number of genes upregulated in white mushroom. (c) Simulation diagram of riboflavin metabolic pathway and riboflavin transport in cells. Arrows represent the direction of the composite, “+” represents that PUT3 positively regulates FlMCH5, and “–” represents that PUT3 negatively regulates PUT2.
Jof 08 00314 g002
Jof 08 00314 g003 550
Figure 3. Functional analysis of FlMCH5. (a) Phylogenetic tree analysis of FlMCH5. (b) FlMCH5 structural–functional domain comparison, where gray represents MFS superfamily. (c) The expression of GUS protein negative control and FlMCH5-GUS in tobacco cells; the cell size in the picture is 20 μm. (d) Riboflavin was extracted from transgenic tobacco leaves; the content of riboflavin was determined to be 3 g in leaves.
Figure 3. Functional analysis of FlMCH5. (a) Phylogenetic tree analysis of FlMCH5. (b) FlMCH5 structural–functional domain comparison, where gray represents MFS superfamily. (c) The expression of GUS protein negative control and FlMCH5-GUS in tobacco cells; the cell size in the picture is 20 μm. (d) Riboflavin was extracted from transgenic tobacco leaves; the content of riboflavin was determined to be 3 g in leaves.
Jof 08 00314 g003
Table
Table 1. Summary of differently accumulation chemical compounds in yellow and white mushrooms.
Table 1. Summary of differently accumulation chemical compounds in yellow and white mushrooms.
ClassTotalDown (in Yellow)Up (in Yellow)Number DifferentRatio (Yellow/White)
Amino acid and derivatives982223450.466
Phenylpropanoids35137200.278
Flavone351898.758
Flavonol162133.691
Flavonoid161120.339
Flavanone112131.392
Isoflavone21122.149
Alcohols155050.005
Polyphenol121120.046
Phenolamides122461.969
Nucleotide and derivates592013331.167
Others2766120.387
Alkaloids288080.103
Carbohydrates2228102.446
Terpene110332.425
Vitamins and derivatives171673.362
Indole derivatives73030.315
Organic acids and derivatives811211231.132
Lipids611218302.274
Anthocyanins21010
Quinones1000/
Sterides5000/
sum573115112227
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gan, X.; Bao, X.; Liu, B.; Li, Y.; Cao, D.; Zhang, H.; Zong, Y. Chemical Constituents and Molecular Mechanism of the Yellow Phenotype of Yellow Mushroom (Floccularia luteovirens). J. Fungi 2022, 8, 314. https://doi.org/10.3390/jof8030314

AMA Style

Gan X, Bao X, Liu B, Li Y, Cao D, Zhang H, Zong Y. Chemical Constituents and Molecular Mechanism of the Yellow Phenotype of Yellow Mushroom (Floccularia luteovirens). Journal of Fungi. 2022; 8(3):314. https://doi.org/10.3390/jof8030314

Chicago/Turabian Style

Gan, Xiaolong, Xuemei Bao, Baolong Liu, Yun Li, Dong Cao, Hg Zhang, and Yuan Zong. 2022. "Chemical Constituents and Molecular Mechanism of the Yellow Phenotype of Yellow Mushroom (Floccularia luteovirens)" Journal of Fungi 8, no. 3: 314. https://doi.org/10.3390/jof8030314

APA Style

Gan, X., Bao, X., Liu, B., Li, Y., Cao, D., Zhang, H., & Zong, Y. (2022). Chemical Constituents and Molecular Mechanism of the Yellow Phenotype of Yellow Mushroom (Floccularia luteovirens). Journal of Fungi, 8(3), 314. https://doi.org/10.3390/jof8030314

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop